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Two recent publications [V. Yakhot, Phys. Rev. E 63, 
026307, (2001) and R.J. Hill, J. Fluid Mech. 434, 379, 
(2001)] derive, through two different approaches that have 
the Navier-Stokes equations as the common starting point, 
a set of steady-state dynamic equations for structure func- 
tions of arbitrary order in hydrodynamic turbulence. These 
equations are not closed. Yakhot proposed a "mean field the- 
ory" to close the equations for locally isotropic turbulence, 
and obtained scaling exponents of structure functions and an 
expression for the tails of the probability density function of 
transverse velocity increments. At high Reynolds numbers, 
we present some relevant experimental data on pressure and 
dissipation terms that are needed to provide closure, as well as 
on aspects predicted by the theory. Comparison between the 
theory and the data shows varying levels of agreement, and 
reveals gaps inherent to the implementation of the theory. 



I. INTRODUCTION 

It is well known that the Navier-Stokes (NS) equa- 
tions of fluid motion can be written in terms of statis- 
tical quantities such as the moments of turbulent veloc- 
ity at several simultaneous spatial points. This statis- 
tical reformulation of the dynamic equations introduces 
extra variables, giving rise to the familiar 'closure' prob- 
lem in hydrodynamic turbulence. If the interest is in 
the small-scale properties of turbulence, it is more ap- 
propriate to obtain equations for the so-called structure 
functions, which are the moments of velocity increments 
Aur — u{x + r) — u{x) over a separation vector, r. 
Such an equation for the third-order structure functions 
is known from Kolmogorov's pioneering work §. This 
equation is unique because it is exact for homogeneous 
turbulence. For the so-called longitudinal velocity incre- 
ment = u{x + r) — u{x), where u is the velocity 
component along the direction of the separation vector, 
the third-order structure function is given by 



{AuD = --{e)r, 



(1) 



where (e) is the average energy dissipation rate. For the 
third-order structure function of the so-called transverse 
velocity increment Avr ~ v{x + r) — v{x), where v is the 
velocity component transverse to the separation vector, 
we also have the result 



{Av^) = 0. 



(2) 



The inertial range is defined by 77 << r << L, where 
T] is the Kolmogorov scale and L is the large scale of 



turbulence. Equations (1) and (2) are both parts of Kol- 
mogorov's 4/5-ths law. 

Recently, Yakhot |Q and Hill ^ have derived dynami- 
cal equations for high-order structure functions. Yakhot 
first derived an equation for the so-called generating func- 
tion Z from which structure functions of all orders can be 
obtained by simple differentiation. Hill used a more con- 
ventional approach to generate structure function equa- 
tions. The equations are new, and it is therefore useful to 
examine if anything further can be learnt about turbu- 
lence through them. However, unlike Kolmogorov's 4/5- 
ths law, these equations are not closed. Yakhot presents a 
"mean-field approach" to obtain pressure and dissipation 
contributions, thereby closing the equations. From these 
closed equations, it is possible to obtain certain small- 
scale properties such as the probability density function 
(PDF) of transverse velocity differences and the scal- 
ing exponents of structure functions of all orders. Our 
goal here is to discuss these new equations for high-order 
structure functions so as to clarify the closure assump- 
tions, and assess the mean field approach by providing 
experimental comparisons for theoretical predictions. 

Since the equations and the procedure for deriving 
them are not yet familiar, we summarize them in Sec. II 
and, for later use, explicitly write them down for struc- 
ture functions of several orders. Section III introduces 
the experimental background needed for our purposes, 
while Sec. IV examines the approximate balance of the 
equations without closure assumptions — mostly to set 
the stage for further discussions. We summarize the 
mean field theory in Sec. V and present in Sec. VI com- 
parisons of its predictions with experimental data on the 
PDFs of transverse velocity increments and their scaling 
exponents. Section VII deals more explicitly with the 
magnitude of dissipation terms, and our conclusions are 
summarized in Sec. VIII. 



II. THEORETICAL BACKGROUND 

A. Brief review of the relevant equations 

Yakhot [H writes the NS equations in terms of the gen- 
erating function Z =< e^'^'^'' >, where Aur is the vector 
velocity difference between two space points Xi and X2 
separated by the vector distance r. The generating func- 
tion Z is constructed, in the spirit of field theories, so 
that its Laplace transform gives the probability density 
function of the velocity differences, and obeys the equa- 
tion 
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dX^drf, 



= If + Ip 



D. 



(3) 



Here If, Ip and D are the (known) forcing, pressure and 
dissipative terms respectively, and are given by 



(4) 



The forcing term in the incrtial range is small and 
may be neglected. The closure of the equation re- 
quires a knowledge of Ip and D. In particular, the 
advection terms are treated exactly. The pressure 
term Ip contains correlation functions of the form < 
{/!^Ur)i{/^Ur)j...{I^Ur)m^(yp) >, and SO One requires 
only the knowledge of the correlation of pressure gradi- 
ent increments and multipoint velocity increments. The 
multipoint energy dissipation function D has a structure 
that depends on details such as the order of the moment 
considered. As we shall see later, its structure resem- 
bles the well-known refined similarity hypotheses of Kol- 
mogorov In any case, the terms needed for the closure 
of Eq. (||) are, in principle, well-defined. 

In homogeneous, isotropic turbulence, the following 
transformation of variables is justified: hi = r, the sepa- 
ration in the direction parallel to the separation vector; 
h2 = the component of A along r; = y/X'^ — 
the component of A in the direction perpendicular the 
separation vector. In these new variables, Eq. ^ for the 
generating function becomes 



dtZ + [did2 + ^-^92 + —d2d^ + 
r r 

(2 - d)h2 ^ ^2 2^r, r , T , r, 
T 03 OslZ = If + Ip + D, 



(5) 



where di denotes the partial derivative with respect to 
hi. In the new variables, the generating function can be 
written as 



Z =< e 



> 



(6) 



where Aur and Avr are the familiar velocity differences in 
the longitudinal and transverse directions, respectively. 
The structure functions are then generated by successive 
differentiation of Z as 

Sn^,n = {{AUrriAVrr) = 8^ Z {r)\,,,^h,=o ■ (7) 

Let us first focus attention on even-order structure func- 
tions for which the dissipation terms are negligible in the 
inertial range. Multiply Eq. (^) by /13, perform SaC^j""^ 
of the resulting equation, and take the limit h2 = hs —> 0. 
We then have 



dS2 



dr 



1 (d-l)(2n-l) _ 

— J2n,0 'J2n-2,2 — \o) 



-(2n-l)(Ap,(A7.)2«-2) 

+ P{1 - COSr/lf)an^dS2n-3,0, 



where d denotes the space dimension, P is the mean en- 
ergy "pumping" rate. If is the forcing scale and a„^d = 

• ^"-"^ n = 1, we obtain the well-known 
relationship between the second-order longitudinal and 
transverse structure functions, as 



55*2.0 d - 1 _ d - 1 

— 7. 1 J2,0 — <-'0.2- 

or r r 

For n = 2, Eq. (^) yields a new relation 
dSifi d-1^ 3{d-l 



(9) 



r, H <5'4.o 

or r 



52,2-3((Apx)(Au)^), 



(10) 



which is exact in the case of incompressible, isotropic 
turbulence in the limit of zero viscosity. To extract fur- 
ther information from this equation, one has to invoke 
some closure. From now on, we consider only three di- 
mensional turbulence unless specified otherwise. If the 
pressure term is small in the inertial range (a considera- 
tion to which we will return momentarily), we obtain 



— 7: 1 '-'4,0 'J2,2- 

Or r r 



(11) 



Following this same procedure, it is now easy to write 
down a sixth order equation with n 3 in Eq. (^). Ne- 
glecting the pressure term again, we obtain 



6,0 



dr 



2o ^lOc 

-J6.O ^4,2- 

r r 



(12) 



In a similar manner, we can extract two additional re- 
lations for fourth and sixth orders. Their corresponding 
approximate forms (again neglecting pressure contribu- 
tions) are 



3^2,2 ^4 ^4 

— 7: 1 — ^2,2 ~ '-'0,4 

Or r 6r 



2,4 



dr 



6 

-i>2A 

r 



6 

or 



(13) 



(14) 



Equations (|Tl|)-(|lJ) were also derived by Hill Q by 
making a convenient change of variables so that the 
scale-dependent parts of the equation can be separated 
from the position-dependent parts and by eliminating the 
position-dependence in the presence of homogeneity. For 
the particular case of isotropy, Hill developed a matrix 
algorithm for solving the system of equations that deter- 
mined all components of the tensorial structure function 
of a given order. The dynamical equations obtained for 
homogeneous and isotropic turbulence confirm Yakhot's 
results. Hill's formulation provides an additional use- 
ful fact that there are exactly two equations relating 
fourth-order structure functions, namely Eqs. ( pi] ) and 
(p^). For the sixth order, his procedure shows that there 
are three equations relating the non-zero components of 
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structure functions, and the third equation is easily gen- 
erated. Again without the pressure terms, this remaining 
sixth-order equation is 



as. 



dr 



-S4, 



-S2 



(15) 



Equations (|llj)-(|_5|) are relationships among different 
components of structure function tensor of the same or- 
der. Without pressure terms, they form a closed set for 
each even order. 

Hill's procedure is convenient for writing down the 
equations for odd-order structure functions. The equa- 
tion for the third-order is the well-known Kolmogorov 
equation (see Eqs. (1) and (2)), which need not be written 
down again. There are three equations for the fifth order, 
which, again without pressure and dissipation terms, are 



dSi 



5,0 



dr 



3,2 



dr 



'5,0 



'3,2 



-S: 



3,2, 



3r 



5*1 , 



5^1,4 _^ 6 

— n ' "-"l,* 

or r 



(16) 



(17) 



(18) 



The « symbol in these equations has to be treated with 
greater caution than for the even-order case because the 
equations for odd-order structure functions contain dis- 
sipation terms as well, see [Q. Therefore the approxima- 
tion imphed in the equations (p^)-(p^) implies the neglect 
of dissipation terms as well. We will assess this additional 
aspect in Sec. VII, but can see by inspection that the ap- 
proximations cannot be correct for at least the last equa- 
tion of the above set. It contains only one component of 
the fifth order structure function, S'1,4, which, when es- 
timated using Kolmogorov's K41 scaling shows that 
the the approximation Eq.(^ cannot hold true. 

A further discussion of these equations, and the of de- 
gree to which they may be reasonable, requires some con- 
tact with experiments. It is therefore necessary to intro- 
duce some basic experimental details at this stage before 
resuming the discussion of the theory. 



III. EXPERIMENTAL CONDITIONS 

The velocity data were acquired by means of a x- 
wire probe mounted at a height of about 35 m above 
the ground on a meteorological tower at the Brookhaven 
National Laboratory. The hot wires were about 0.7 mm 
in length at 5 /zm in diameter. They were calibrated just 
prior to being mounted on the tower, and operated on 
DISA 55M01 constant-temperature anemometers. The 
frequency response of the hot wires was typically good 
up to 20 kHz. The voltages from the anemometers were 
low-pass filtered and digitized. The low-pass cutoff was 



never more than half the sampling frequency f^. The 
sampling rate was adequate to resolve most of the scales, 
including dissipative ones. The voltages were converted 
to velocities in a standard way through the standard non- 
linear calibration procedure. The mean wind velocities, 
roughly constant over the duration of a given data set, 
ranged between 5 and 10 ms~^ in the entire series. The 
usual procedure of surrogating time for space ( "Taylor's 
hypothesis" ) was used to obtain the mean dissipation rate 
(e) and an estimate for the Kolmogorov scale rj. 

The relevant details of the data analyzed here are given 
in [^[0]. Briefly, the Taylor microscale Reynolds number 
is 10,680, the large scale L is about 42 m, Kolmogorov 
scale 77 is 0.44 mm, the scaling range according to the 
linear part in the third-order structure function is con- 
servatively between 0.01 m and 0.2 m (although it could 
be stretched in either direction by factors of 2). This will 
be regarded as the operational definition of the inertial 
range. 

The theoretical development discussed here is meant 
for locally isotropic turbulence appropriate to asymptot- 
ically large Reynolds numbers. The present measure- 
ments are indeed at large enough Reynolds numbers, but 
it is not obvious that the small scales are isotropic. In- 
deed, we have used similar data before js]]^ to extract 
anisotropic parts of the structure functions by perform- 
ing the S0(3) decomposition. Thus, a few explanatory 
words regarding the degree of isotropy are appropriate 
here. 

Figure |^ shows the longitudinal spectral density En 
and the transverse spectral density E22 for one of the 
data sets. If strict isotropy prevails in the inertial range, 
the ratio E22/E11 should be 4/3. The data show that 
the ratio, while varying slowly in the inertial range from 
a somewhat smaller value than 4/3 to a somewhat larger 
value, is not far from being 4/3. Similarly, for third-order 
structure functions, we have the exact result jl0| that 
the ratio S'i,2/<S'3.o should be 1/3. Measurements show 
that this ratio, which is indeed reasonably constant in the 
inertial range, has a magnitude of about 0.42 (see Fig. ||). 
There is undoubtedly some degree of anisotropy in the 
inertial range, and so the conclusions are to some degree 
affected by this artifact. To pursue this issue further, 
we note that 83^0 has the expected value of 4/5, and 
So,3, which should be zero exactly for the isotropic case, 
is of the order of 0.1 83^. The anisotropy is thus not 
large, and the transverse component is perhaps the more 
anisotropic. This conclusion is consistent with [p|JTT|. 

In summary, then, we regard for present purposes the 
departures from isotropy to be relatively small in the in- 
ertial range, especially for even orders, and expect that 
the results will not be qualitatively affected by their pres- 
ence. The quantitative effect is not easy to ascertain, 
though it is likely to be modest from the considerations 
just cited. 
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IV. THE BALANCE OF THE APPROXIMATE 
EQUATIONS 

Using the experimental data just described, we calcu- 
late the left-hand side (LHS) and right-hand side (RHS) 
of each approximate equation of Sec. II, and obtain the 
relative size of the difference (LHS - RIIS)/LHS; for even 
orders, this difference is an estimate of the magnitude of 
the neglected pressure terms. Figures H, l 4[ ||, M and mdis - 
play the results for Eqs. (|ll|), (|l|), (|l3^ ^ll) a^d^ll) , 
respectively. The absolute value of the relative difference 
is also shown in each case. In the inertial range, all four 
equations seem to balance reasonably well without the 
pressure terms, their contribution being of the order of 
10% of the LHS for all equations, growing larger as r 
becomes larger as well as smaller. To that extent, the 
closure of Eqs. (pl])-(p^) in the inertial range seems to be 
approximately justified. 

The situation with respect to odd-order structure func- 
tions is somewhat different. Odd-order moments do not 
converge as well as even-order moments do, and so we 
may expect significantly more scatter in experimental 
plots. The terms for the three equations (16)-(|8|) are 
shown in Figs. ^10. Despite the large scatter, it is rea- 
sonably clear that the relative difference of LHS and RHS 
in Fig. I is of the order of 20 — 30% in the inertial range. 
In Fig. |, the relative difference is perhaps larger, up 
to 50% at places. For Eq. (p^), as was expected and 
remarked upon earlier, there is no balance at all (see 
Fig.0). 

Before examining these issues, it is worth commenting 
on the intriguing observation that the pressure effects 
are small in the inertial range at least for even- order 
structure functions. This feature suggests that the op- 
erating physics there might have some relation to that 
of the forced Burgers equation. It further suggests that 
the pressure terms may become effective only in the dis- 
sipation range (and also for large scales, which we will 
not consider here). It is reasonable to suppose that vor- 
tex structures (qualitatively like wing-tip vortices) are 
generated as soon as the pressure terms are activated. 
Perhaps the small-scale vortex tubes (also called worms, 
see [^2[|3)), are a result of this effect. Clearly, the va- 
lidity of the pressureless physics is limited because the 
scaling exponents saturate for the forced Burgers equa- 
tion, whereas there is no evidence that the saturation 
occurs for three-dimensional turbulence (at least for mo- 
ments up to order ten; see Sec. VI). In any case, the 
smallness of pressure terms suggests that the intercom- 
ponent energy transfer (for which they are responsible) 
is small, and that any remnant anisotropics at the large- 
scale end of the inertial range tend to persist, or at least 
not diminish rapidly, as the scale size decreases through 
the inertial range. There is growing evidence from other 
measurements that this might indeed be so |9|,[T^. 



V. A MEAN FIELD THEORY 

The previous section has shown that the imbalance in 
the approximate equations for even-order structure func- 
tions, caused entirely by the neglect of pressure terms, 
is of the order 10%. The imbalance is larger for odd- 
orders, for which it is to be remembered that dissipation 
terms are also important. Thus, ignoring pressure and 
dissipation terms is not an option in general, and it is 
clear that one must make a plausible theory for them. 
This has been attempted in the mean field theory of [|. 
Though our interest and contributions are primarily in 
the experimental assessment of the theory, it is helpful 
to provide here a summary of the theory itself — if only 
to clarify the motivation for the experimental tests per- 
formed. We shall focus on the essence of the physical 
arguments, rather than on analytical details. 



A. General remarks 

Mean field theories provide approximate means of de- 
scribing a thermodynamic system by supposing that each 
'particle' in a many-body system moves in the 'mean' 
field of all other particles in the system. This is op- 
posite to the situation in which only nearest neighbor 
interactions matter. More formally, attribute to the sys- 
tem an order parameter that is zero when the sys- 
tem is ordered and becomes increasingly non-zero with 
increasing disorder. If the fluctuations in the order pa- 
rameter are small, then it may be replaced by a spatially 
uniform average value. The mean field approximation 
implies infinite range interactions; while this cannot be 
realized in practice, the order parameter in many ther- 
modynamic systems could become arbitrarily small as 
the temperature approaches a phase transition value, Tc. 
The Ginzburg-Landau theory makes use of this feature 
to propose a description of the free energy and to derive 
critical exponents at phase transitions. In general terms, 
the free energy F{(j3, T) is expanded in powers of (j) as 



F = FQ + A(t)- 



(19) 



where A, B and C are functions of T. Near the criti- 
cal point in the T-space, where T Tc, the expansion 
can be truncated at the lowest order terms in (jj. The 
expansion then provides a qualitative description of the 
thermodynamic processes; in practice, this mean field ap- 
proach may work even far from the critical point. 

Strictly speaking, a mean field theory may not ap- 
ply to turbulence where quantities such as the free en- 
ergy and order parameter cannot be defined unambigu- 
ously. In Yakhot's theory, the idea is carried over quali- 
tatively by identifying a small parameter in some regime 
and expanding other dependent quantities around that 
small parameter. The 'phase transition' considered is 
the change of sign in energy flux that occurs in going 
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from two-dimensional (2d) to three-dimensional (3d) tur- 
bulence. It is understood from Kolmogorov's equation 
for the third-order structure function that the energy 
transfer is from the small to the large scale in 2d turbu- 
lence, and vice versa in id turbulence. It is assumed that 
this change is continuous, changing sign at some critical 
dimension dc — analogous to the critical temperature Tc 
in thermodynamic phase transitions. In 2c? turbulence, 
the dissipation is negligible for high Reynolds numbers 
(because the energy ultimately concentrates in the large 
scale). In 3d turbulence, on the other hand, the dissi- 
pation is the key to energy transfer from large to small 
scales. The hope, then, is that both pressure and dissi- 
pation can be expanded in terms of a small parameter in 
the vicinity of dc- 



from Eqs. (13)-(15). This is in excellent agreement with 
the results of numerical simulations of jl^. Thus, we 
might conclude that a plausible mean field expression for 
the pressure contribution exists for 2d. 

The next crucial assumption of the theory is that the 
above form of the mean field approximation is applica- 
ble also for 3d turbulence. The rationale is not easy to 
articulate, especially because, unlike in 2d turbulence, 
the PDFs of /S.Vr possess stretched-exponential tails in 
3d turbulence . We shall provide some statements of 
mild justification subsequently, but emphasize that the 
validity of this assumption has to rest on the basis of the 
agreement, or lack thereof, with experiments. 

C. The small parameter and the dissipation term 



B. The pressure terms 

In turbulence dynamics, transverse structure functions 
do not participate in energy transfer. Yakhot therefore 
regards the fluctuations in the transverse velocity incre- 
ment Ai;^ as small — in effect, if not in actual fact. It is 
known from numerical simulations [ p^ as well as exper- 
iments of the inverse cascade in 2d turbulence that 
is almost exactly gaussian. The absence of inter- 
mittency makes it plausible to regard the fluctuations as 
"small" . We shall therefore consider the 2d case briefly. 

The key step for further analysis is the introduction of 
a conditional expectation of the pressure gradient incre- 
ment for a fixed value of Au^, Aw^ and r as 

{dyp{x + r) - dyp{x)\/^Ur,/^Vr,r) (20) 
«^«,n.n(r)(Au,)"(A«,)". 

This is related to the needed correlations in Ip which is 
of the form 

( (dyp{x + r)- dyp{x)){^UrY{^VrY) (21) 

{dyp{x + r) - dyp{x)\/^Ur, Avr,r)AuPAv^ 

X P{Aur, Avr,r)d{Aur)d{Avr)d^r. 

The use of the conditional expectation provides a tool 
for expanding the pressure terms in terms of the "small 
quantity" Au,.. Now, in the spirit of the Ginzburg- 
Landau expansion, only the lowest order terms in Av^ 
are retained (corresponding to Au^Avr and Avr). The 
prefactors of the expansion are constrained by the incom- 
pressibility condition and by the dimensionality of space. 

By substituting in Eqs. (9)-(ll) the pressure term de- 
rived from the conditional expectation value, and assum- 
ing the exponents to be given by from Kolmogorov's K41 
scaling arguments Yakhot concludes that the high- 
order even moments are consistent with gaussianity. The 
argument is circular but internally consistent. The gaus- 
sianity of the transverse increment Avr is then deduced 



We need to consider the dissipation term before turn- 
ing to experiments again. In the inverse cascade range in 
2d turbulence the dissipation term D can be set to zero 
because the flow evolution is towards larger and larger 
scales. However, D is central in 3d turbulence, and it 
is known that dissipation fluctuations are immense at 
high Reynolds numbers [0. The objective in a mean 
field approach is to locally smooth out the fluctuations, 
through some procedure such as Obukhov's [|l8|. For 
closure, there is a need to relate this coarse-grained dis- 
sipation field to velocity fluctuations, analogous to that 
employed in Kolomogorov's refined similarity hypotheses 
1^,^. Yakhot 's theory is similar in spirit but the details 
are different, as we shall illustrate. 

Let us denote a coarse-grained velocity field for a given 
spatial scale r by This will be assumed to be the 
same as Avr. Certain one- loop calculations of Yakhot 



and Orszag |21| give the effective viscosity as 



(d- 



)iN{ey)i (22) 
higher order nonlinear terms 



where is a constant that depends weakly on the space 
dimensionality, d, and Sr is the dissipation rate coarse- 
grained on the scale r. If we ignore nonlinear terms, this 
equation provides a natural definition of r^, the charac- 
teristic time for the field Vr- 

There is no obvious justification for ignoring the higher 
order nonlinear terms which, in 3d turbulence, are typi- 
cally 0(1), nor in assuming that is small compared to 
Vr/r. However, if we assume that the theory can be an- 
alytically continued into non-integer dimensions between 
2 and 3, a suitable small parameter can be generated as 
follows. The time scale characterizing the interaction of 
a scale r with all other scales less than r is the so-called 
eddy turn-over time, or the time taken for energy trans- 
fer to occur between r and the Kolmogorov scale rj. One 
may use K41 to estimate this time scale. The process 
of energy transfer can be thought to consist of two dis- 
tinct steps, one involving nonlinear transfer across scales 
without any pressure effects, and another involving the 
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relaxation due to pressure effects. In 3d turbulence, these 
two steps are parts of the same simultaneous process, 
so the time scales associated with them cannot be sep- 
arated. But, if, as one approaches dc, it is increasingly 
true that the pressure effects are small except when scales 
of the order 77 are reached, the two time scales involved 
could become disparate, and the relaxation due to pres- 
sure terms enters the picture only at the smallest scale 
and can therefore be assumed to be fast. Then the di- 
mensionless ratio qr = Tr/9r, where is the time scale 
for relaxation effects and Or is the time scale for energy 
transfer, would be a small parameter. 

Using this basic idea and his one- loop calculations , 
Yakhot deduces the following results: 

4 = 2.56, qr^{d- Vr^{d~ 

and so forth. The notion of a critical dimension is not new 
(see iQ), though the estimates for it in |H and [|2| are 
substantially different. The precise numbers and powers 
in the above equation depend in detail on the approxi- 
mation made to compute them, and are presumably not 
final; they cannot, in any case, be verified experimentally 
near the critical dimension. Here, we merely wish to draw 
upon the general idea of a critical dimension near which 
a small parameter can be defined, and in whose vicinity 
the energy piles up (as shown by the last of the three 
relations above: the energy is being pumped at a con- 
stant rate but is being transferred neither upscale as in 
2d nor downscale as in 3d). These ideas allow Yakhot 
to truncate the effective viscosity and write the dissipa- 
tion Er in terms of the lowest order terms in terms of the 
coarse-grained velocity field in the vicinity of dc'. 

Sr ~ -l-^[V„V^jil + Old ~ d,))] (23) 

Perhaps two additional remarks might be usefully 
made. First, the coarse-grained velocity fluctuations be- 
come very large as the critical dimension is approached, 
yet it may seem that the mean field approximation pro- 
posed for pressure terms assumes that fiuctuations are 
small. To avoid confusion, it is important to keep in 
mind the distinction between fluctuations in longitudinal 
and transverse velocity increments. The velocity scale 
that blows up is related to energy transfer, and hence 
the longitudinal velocity component, but the component 
whose fluctuations are supposed to remain small is the 
transverse velocity. The sense in which those fluctuations 
are small is unclear (because they too are intermittent in 
3d, see fl^), but the fact remains that it takes no part 
in energy transfer and so its effects are thought to be 
"small" in some rough sense. Since the pressure effects 
are small, the intercomponent energy transfer is inhib- 
ited, and so, once fluctuations in Avr are small at some 
scale, they will presumably remain small at others as well. 
Secondly, in order to be able to truncate the energy dis- 
sipation, the higher order viscosity terms have to decay 



faster than the rate of blow-up of velocity fluctuations. 
This is indeed the case above. 

Now, keeping in the mind the symmetries of the NS 
equations, the simplest form for the Vr contributions to 
the dissipation rate is 



(e) w c(d)AurAwi 



dAu, 
dr 



(24) 



The coefficient c(d) must refiect the change in going from 
2d to 3d (zero dissipation to finite dissipation) . This may 
not be a smooth change (as in second-order phase tran- 
sitions) because c(dc) could well be singular (as in first- 
order phase transitions). Yakhot assumes, however, that 
it is 0(1) for d — dc > 0. Then D takes on a form similar 
to Kolmogorov's refined similarity hypothesis, relating 
{Sr) with the third-order longitudinal structure function: 



D w c{d)h3dh2dhadrZ 



(25) 



c{d)hl{AurAvr 



^^^^h2Aur.+h3Av,. 

dr 



+ terms neglected. 

This enables the closure to be complete. 

A non-trivial difficulty is the testing of the theory in 
non-integer dimensions near dc- At present, the conse- 
quences of the theory can only be tested in 2d or 3d. The 
extrapolation to non-integer dimensions is not an intrin- 
sic limitation of the theory, but reflects the lack of exper- 
imental ingenuity at present. It must, however, be noted 
that in shell models where an interaction parameter can 
be tuned to change the direction of energy transfer, one 
can make more reasonable contact with the theory. Such 
comparisons have been attempted recently [^4| and the 
results are encouraging. We have already noted that the 
conclusions of the theory are consistent with experiments 
and simulations of 2d turbulence. Wc shall examine in 
the rest of the paper the extent to which the predictions 
of the theory are applicable also to 3d. 



VI. COMPARISON OF THE THEORY WITH 
MEASUREMENTS IN THREE-DIMENSIONAL 
TURBULENCE 

A. Probability density function of transverse 
velocity increments 

When the forms of Ip and D from previous sections 
are substituted into the full structure function equations, 
one can generate the following equation for the PDF of 
transverse velocity increments P{Avr,r) = P{V,r): 



dP 

dr 



1 + 3(3 d 
3r dV 



VP 



(26) 



Here /3 oc c(d) . This equation is linear and can be solved 
in principle, but the solution has no simple analytic form. 
(For some discussion of this aspect, see [H.) For small 
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V, however, the equation admits a solution of the type 
r''F(T//r«) with 



P{V = 0,r) a 

where, from Yakhot's theory, 

^ 1 + 3/3 
"^"3(1-/3) 



0.4 



(27) 



(28) 



and f3 ~ 0.05. Using the fact that the variance ay of V 
is expected to vary as r" '^^ (see 0), we have the result 



P{— 0,r) = ayPiV 

(JV 



0,r) 



(29) 



We shall now test these predictions. 

The precise measurement of the peak value of the 
PDFs from the data must be done carefully because it 
is sensitive to the bin width chosen around ^ = 0. In 
our measurements, the bin width around = was 
gradually refined until the PDF value at the origin no 
longer depended on the bin size. Figure 11 shows that 
P{y ^ 0) at r « 100?7 asymptotes to a value of 0.64. 
The sharp ascent of the numbers for very small values 
of the bin width is an artifact of the extreme narrow- 
ness of the bin width, which results in false values when 
normalizing. The procedure was repeated for several val- 
ues of r. Figure H shows the properly normalized PDF 
values for V = Q for different scales r ranging from the 
Kolmogorov scale -q to the large scale L. The scaling ex- 
ponent for this quantity is k, —0.065 in the inertial range, 
numerically about 25% larger than the theoretical value 
of -0.05. With this experimentally derived scaling expo- 
nent, we can evaluate that (3 « 0.058 compared to the 
theoretical value of 0.05, a 16% difference. Figure |l3| also 
shows that the form r'^PiV/r'^) is essentially constant for 
small V . 

One can obtain the form of F for large V hy a, steepest 
descent approximation (see |2^ for more details). The 
result from iHl is 



F cx 



1 



: exp 



Ml: 

n{r) 



(30) 



]/2,7?^/r«(i-2/3)/i-/3^ and n{r) = 
(1-/3)2 1 -^-^^v I (The corresponding expressions in ||] 
are printed incorrectly.) The prefactor of Eq. ( pO| ) is pos- 
sibly r-dependent. Equation (^0|) can be re- written as 



where ^ = 
^^|ln(r/L) 



(-rj(r)ln(P(V^,r)r''VfX^)°'^ « ln(0- 



(31) 



Figure |TJ shows plausible linear behaviors for the tails of 
the PDF in the proposed logarithmic units of Eq. (31). 
There is, however, evidently still some r-dependence that 
precludes their collapse. We recall that corrections to 
steepest descent approximations are often logarithmic, 
but are difhcult to calculate here analytically. We assume 



a dependence of the form (ln(r/L))'>' for the proportion- 
ality factor. Figure |l^ shows a replot of the data with 
the additional factor of (ln(r/L))^ multiplying the PDF. 
The exponent 2 was chosen because it collapses the data 
best in the inertial range. (The one separation distance 
that does not collapse belongs to the dissipation range.) 

Our main conclusion so far is that the mean field mod- 
els for pressure and dissipation terms provide a way for 
closing the PDF equation, and for solving it for the lim- 
iting situations. The prediction is that, to first approxi- 
mation, the tails of the PDF of are lognormal. The 
experimental data suggest that this might be so, but that 
an r-dependent contribution is missing. It is at present 
not clear whether this missing aspect is merely a correc- 
tion to asymptotics, or corresponds to additional terms 
in the mean field expansion, or is even more fundamental. 



The scaling exponents and the prospect of their 
saturation 



Seeking the solution to Eq. (g6|) under the K41 con- 
straint for the third-order structure functions and assum- 
ing 5*0, n oc r'-" , Yakhot obtained the following formula for 
the structure function exponents: 



Cn 



n(l + 3/3) 
3(1 + (3n) ■ 



(32) 



Table I and Fig. |l^ show the calculated exponents and 
compare them with those obtained from the Direct Nu- 
merical Simulations, or DNS, data as well as exper- 
iments. The agreement is good for all orders, perhaps 
slightly better for the DNS data for high order exponents. 

Using probability density functions to define the sta- 
tistical quantities, we have (putting U — Au^) the con- 
ditional expectation value of for a fixed value of U , 
Q2{U), as 



^2n-2,2 = / P{U)U'''-'Q2{U)dU. 



(33) 



See Sec. IV of [0. The Kolmogorov scaling will hold 
(by dimensional arguments) for Q2{U) oc U^. On the 
other hand, saturation of exponents, C2n —> constant as 
n — > oo, is possible for Q2{U) oc for 5 <2 and U large. 
We present the conditional statistical quantity Q2{U) as 
a function of r in Fig. ^ It is not clear if the trend for 
large U is in agreement with the saturation condition. 
There is a very small range of U towards the tails which 
seems to vary as t/^/ln(C/^) but this is not conclusive. 
There might also be the influence of anisotropy in the 
PDFs, as is evident, for example, in the asymmetry of 
the joint PDFs, which in turn could change the nature 
of the tails of the conditional statistics. 
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VII. REMARKS ON THE MAGNITUDE OF 
DISSIPATION TERMS 

It is helpful to recall that the full form of the equa- 
tion relating even order transverse moments to mixed 
moments of the same order in 3d is 



2,2n 



dr 



2n _ 2 

i>2,2n — 



2n S, 



0,2n+2 



2n+l r 
- 2n{Py,Au{Avf''-^), (34) 



where Vyv = dyp{x + r) — dyp{x). The subscript y de- 
notes the component of the pressure gradient while the 
subscript v reminds us that this form of the pressure may 
be expanded in Av^. In the inertial range, the dissipa- 
tion term is small (because the order of the structure 
functions is even) and the forcing term negligible. 

Compare this equation to the one relating the odd or- 
der mixed moments 



dSi,2n ^ 2 + 2n 
dr r 



Si^2n = -2n{Vyy{Ai 



')-D. (35) 



The dissipation term is present in this case, and the pres- 
sure term shares the same derivative factor Vyv found in 



Eq. (34). Keeping the first two terms in the mean field 
expansion (see the first equation in Sec. VB), we have for 
that term 



AUr-AVr AVr 

B- 



[Pr) 



(36) 



where H and B are unknown constants, and P, as be- 
fore, is the rate of forcing. The second term is chosen 
in order that the pressure term in the third-order equa- 
tion VyvAvr = 0, which is a K41 constraint. The two 
constants H and B are then related through 



HSi,2 — BSq,; 



^1/3 



(37) 



If we insert the model for the pressure into Eq. (^4|) for 
n = 1 (say), we can extract the value of the constant 
H from experimental data on the left and right sides of 
that equation. With this value of H in Eq. ( p5| ) we have a 
complete expression for the pressure term. The remain- 
ing imbalance, if any, must be the dissipation term. 

To evaluate the magnitude of the dissipation terms, let 
us return to Eq. (|34|), for which, with 



B- 



AVr 



(Pr) 



(38) 



we have the following result: 



2n{VyvAur{Avr)^"-^) = -2n{-H 



S2,2n „ Sl^2n 



B- 



{Pr)i 



„ rj,S2,2n 1 5'l,2'S'i 2n s .„„n 

2nH{ . 39 

r r io,2 



The first equality in En. (^^ uses Eq. (^6|) and the sec- 
ond equality uses Eq. (^7|) . We extract the value of H by 
substituting Eq. (^) in Eq. (^J) and solving for it from 
experimental data for the case n = 1. Figure |l^ shows 
H computed in this manner. The value ranges between 
0.3 and 0.1 in the inertial range with significant scatter 
for larger scales. The value of H can also be cornputed 
from the sixth-order equation (i.e., n = 2 in Eq. (|34|)) in 
this same way. Though the statistics are not as nicely 
convergent, the result is that H ss 0.28 with deviations 
of the order of 0.1. As a further check, we can compute 
H from the above equation assuming the scaling expo- 
nents ^2,2 = Co, 4 — 1-26 and taking from the data the 
ratio S'o,4/iS'2,2 ~ 4.5 in the inertial range. We obtain 
H « 0.37. Thus, the precise value of H appears to de- 
pend on the order of the moment and on the scale range 
but, since the pressure term is relatively small in the in- 
ertial range, the exact choice of H may not be critical in 
determining the dominant dissipation term. In any case, 
the uncertainty in these estimates does not allow us to be 
too definitive, and so we shall proceed with an average 
value of 0.25. 

Our proposal is to substitute this value of H in Eq. ( |35| ) 
for n = 2. In detail, that equation is 



2n{Pyy{Avrr''-'-) = -2n{-H 



Si.:. 



B 



So,: 



{Pr)i 



2nH{ 



Si,: 



1 '5'l,2'S'o.; 



Sq,: 



(40) 



This pressure term is seen to account for only about 10% 
of the imbalance of Eq. (35), as shown in Fig . 09^ . We now 
substitute the pressure contribution of Eq. (|40|) back into 
Eq. ( pS) in order to estimate the only unknown term D. 
Figure 2^ shows that the dissipation term so extracted 
dominates RHS, and is much larger than the pressure 
term. The dissipation term alone is comparable to the 
entire LHS of the equation. 

There is another equation involving fifth-order struc- 
ture functions that contains the same form for the pres- 
sure gradient as Eqs. (|4|) and (^). It is Eq. (|l^) in- 
cluding the pressure and dissipation terms which, in full, 
reads as 



a c A o 

' -S3,2^^S,,4-2{Py,{AUrfAVr)+D. (41) 



dr 



We now follow a similar procedure as before. For 
Eq. (41), the pressure term according to the mean field 
model is 

- 2(n.(A«.)^Azv) = -2i-H^ B^^) 

= 2ff(^-i%^). (42) 
r r 5o,2 

Figure ^ shows that the RHS of Eq. ( p^ ) balances 
the LHS up to about 80% in the inertial range. The 
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imbalance is due to a possible mix of pressure and dis- 
sipation. The pressure term computed from Eq. ( ^2|) is 
shown in Fig. ^ It makes a 10% contribution in the in- 
ertial range. The dissipation term, being the remainder 
{D ^ LHS - RHS 



23 



while it 



Ip), is plotted in Fig. 
shows significant scatter, it is clearly small in the inertial 
range (of the order of 15% or less) while increasing, as it 
must, toward the dissipative scales. 

From the above two examples it appears that the over- 
all order of the structure function (in this case the fifth) 
is not enough to prescribe the relative importance of the 
pressure and dissipation terms. The equations that relate 
diff'erent components of the fifth-order structure function 
tensor: Eq. (^) shows different ratios of pressure and dis- 
sipation terms from Eq. (|3^), whereas the overall order 
of the structure function in both cases is 5. The former 
equation seems to balance more or less without pressure 
and dissipation while for the latter equation the pressure 
term and, particularly, the dissipation term prove to be 
essential. 

There is one further detail that needs to be mentioned 
for completeness. This concerns the relative magnitudes 
of the H and B terms in Eq. (|4^). The ratio of the 
B term to H term in Eq. (42) is about 10%, while this 
ratio is of the order of 1% for Eq. (|^ and approximately 
80% for Eq. (^5|). The conclusion seems to be that the 
relative importance of the i?-term depends on the order 
of the structure function as well as the component of the 
structure function being considered. 



VIII. CONCLUDING REMARKS 

Our experimental results are assessed in the context 
of a mean field model due to Yakhot. The model allows 
us to write the pressure terms which we cannot measure 
directly, in terms of the velocity structure functions that 
we can measure. (The pressure terms appear here in a 
different form from those used in turbulence modeling, 
and so the value of the present work to that endeavor is 
unclear.) Among the assumptions made, the most drastic 
one is the use of the same pressure model for 2c? and 3d 
turbulence. 

Nevertheless, if we adopt the pressure model in 
Eq. ( 34 ) in which the dissipation terms are thought to be 
negligible (see ||^ and for symmetry and asymptotic 
arguments as to why this might be so), the coefficients 
H and B can be obtained, and thus the pressure terms 
can be modeled. We can now proceed to analyze odd- 
order equations that have the same structure for pressure 
terms. Since the pressure term is known, we can deduce 
the only remaining term, namely the dissipation. For 
one equation (p5|), the dissipation term is of the order of 
80% of the balance. Another dynamical equation (^ ) 
for the same order of the structure function has a differ- 
ent structure, and there, the dissipation term is relatively 
small. This is a new and interesting statement about the 



inertial range dynamics, but its validity depends on the 
pressure model used. At least one outcome of the calcula- 
tions is tautologically correct: in all the cases considered 
here, the dissipation range is always dominated by the 
dissipation term D. 

Yakhot 's theory postulates the existence of a critical 
dimension, dc- This, in itself, is not implausible [ p2[ . 
However, the analytic structure of the NS equations in 
the neighborhood of dc and the extent of the neighbor- 
hood remain unclear. The theory yields certain expo- 
nents for the vicinity of dc, but the details on which they 
are based need closer scrutiny; at least to us, some of the 
steps remain unclear. Thus, while the numerical values 
of the exponents, as well as that of dc itself, are proba- 
bly not to be taken literally, we should be interested in 
drawing some qualitative conclusions. 

Such conclusions come from a few independent sources. 
First, the prediction of the theory for the PDF of Avr for 
2d turbulence is in good agreement with simulations and 
experiments ||l^,|l^. Second, the conditional expectation 
of the pressure terms in 3d simulations [ p6| appear to 
follow the mean field theory, at least for modest values 
of the velocity increments. Third, shell model calcula- 
tions ||2^ show that the behavior expected near the crit- 
ical dimension can be observed as one varies a coupling 
parameter. Finally, the present comparisons with exper- 
imental data at high Reynolds numbers reveal that the 
scaling of the PDF of Avr for small and large Avr are 
in some measure of agreement with the theory. All these 
are positive developments. However, since many details 
are unclear, it remains to be seen as to whether the the- 
ory will evolve into a rational framework. For now, we 
find it to be remarkably interesting and worthy of some 
attention. 

We thank V. Yakhot for numerous helpful discussions 
without which this paper would not have been written. 
However, he should not be held responsible for our in- 
terpretations of the theory. We also thank R.J. Hill and 
J. Schumacher for helpful comments, and the latter also 
for reading the manuscript carefully. The work was sup- 
ported by the US Office of Naval Research. 
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TABLE I. Comparison of exponents from the DNS data 
and the experiment (both using ESS) and Yakhot's formula, 
Eq. (H). 
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FIG. 1. Longitudinal and transverse spectral densities, 
Taylor microscale Reynolds number ~ 19, 000. 
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FIG. 3. The terms of Eq. (|ll|) with their relative difference. 
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FIG. 2. The ratio of 51,2 to 53,0 showing some deviations 
from the isotropic value of 1/3 in the inertial range. In this 
figure and elsewhere, the horizontal bar in the figure repre- 
sents a conservative extent of the scaling region determined 
from Kolmogorov's 4/5-ths law. 
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FIG. 4. The terms of Eq. (02) with their relative difference. 
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FIG. 9. The terms of Eq. (0) with their relative difference. ^^G. 11. The value of the PDF as VIgv ^ for r ^ IOOt; 

for different binwidths. 
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FIG. 10. The terms of Eq. (|l|) with their relative differ- 
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FIG. 12. Log-log plot of the peak values of the PDFs of 
transverse velocity increments. The line indicates a slope of 
-0.065. 
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FIG. 13. Log-log plot of the near-peak values of the PDFs 
of transverse velocity increments. The collapse of the data 
occurs for the normalization r'^P{V/r"') where k = 0.065. 
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FIG. 15. The collapsed tails of the PDF from the previous 
plot, taking into account an ad hoc prefactor of (ln(r/I/))^. 
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FIG. 14. The tails of the PDFs in the form required by the 
theory (see Eq. (^l|). Lognormality of the tails requires that 
the data follow straight lines. 
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FIG. 16. Comparison of the DNS and experimental values 
of C^n with those from Eq. (^2|). The numbers in the plot are 
for longitudinal structure functions, but these are identical to 
those for transverse structure functions in the isotropic sector 
see 
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FIG. 17. Q2{U), the conditional expectation value of V'^ 
on U for various different r, on a log-log scale. The upper 
solid line indicates /log{U^), the lower solid line indicates 
the scaling slope. 
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FIG. 19. The LHS and pressure terms of Eq. ( psj ) with 
n = 2 and H — 0.25 in the pressure model; the ratio indicates 
that the pressure term as computed from the model is about 
10% of the balance. 
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FIG. 18. The coefficient H of the pressure term computed 
by substituting Eq. Q in Eq. with n = 1. 



FIG. 20. The terms of Eq. |35| with n = 2, H ^ 0.25; the 
difference LHS-RHS indicates the magnitude of the dissipa- 
tion term. 
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FIG. 21. The LHS and RHS of the approximate equation, 
Eq. (|l7|), without pressure or dissipation terms. The RHS 
balances the equation up to about 80% in the inertia! range, 
less as we approach the dissipation and large scales. 
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FIG. 23. The LHS and dissipation terms D of Eq. 
with the dissipation computed by subtracting the RHS and 
the pressure (see Fig. from LHS; the ratio of D to LHS in- 
dicates a large scatter but shows the dissipation contributes in 
the inertial range only to about 10%, but that it significantly 
increases towards the dissipative scales. 
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FIG. 22. The LHS and pressure terms Ip of Eq. @ with 
their ratio showing that the pressure only accounts for about 
8-10% of the balance, less as one approaches the dissipative 
scales and large scales. 



16 



